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BECHTEL CORPORATION (A) 
Process Analysis and Optimization of a Desalting Plant Design 

Mr, Rufus Crawford, a young chemical engineer, had been hired by the 
Bechtel Corporation in 1964 to assist in process analysis and cost studies 
of large multi-stage flash saline water conversion plants. An earlier 
Bechtel study was sponsored by the Office of Saline Water (OSW 116) to de- 
termine how much costs for an M.S.F. plant could be reduced by increasing 
the size of the operation and using exhaust steam from a gas-fired power- 
generating station to heat the sea water. 

Rufus said that they had tried various optimization techniques using 
programs they had obtained from the National Computer Library. He said, 
"It's hard to pick up the logic of those library computer programs without 
spending a great deal of analysis time/* 

Thus, Crawford was enthusiastic when Dr. Lionel S. Galstaun and 
Mr. Richard L. Clark called him into Dr. Galstaun's office to ask him if 
he would like to work with Dr. Douglas Wilde of Stanford in the development 
of a new optimization program. 

The Water Resources Section of this large engineering and construction 
firm in which Crawford worked had just received management approval for an 
in-house sponsored project to study the feasibility of large scale flash 
distillation in order to (1) improve the company's competence in desalting 
technology and (2) develop experience in computer-aided optimization. 

Although the company was well experienced with desalting on a pilot, 
plant scale, large plants .would be something new. Tradeoffs to be made 
would include such things as energy required versus investment in heat 
transfer area. What temperatures and pressures to operate at, and what 
kinds of water treatment (sulphuric acid or phosphate additives?) to use 
for suppression of scaling were among the questions in need of answers. 

Crawford saw his opportunity to really answer some of those questions 
he had been asking himself about those computer programs and gladly accepted 
the assignment. 
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Prepared at Stanford University during the 1966 National Science Foundation 
Summer Institute conducted by the Design Division of the Mechanical Engineering 
Department, The initial contact with Bechtel was made by Karl H. Vesper . The 
case was prepared by Joseph W. McCutchan, U.C.L.A., Robert P. Vail, University 
of Arizona, and Edward N. Ziegler, Polytechnic Institute of Brooklyn, parti- 
cipants in the Summer Institute. The assistance of Mr. Rufus Crawford was 
indispensable since the case really relates his experience. The cooperation 
of Mr. James Appleton and the permission of the Bechtel Corporation to use 
this case for educational purposes is greatly appreciated. 
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THE EECKTEL CORPORATION (A) 

Applied Technology Group oi the 
Scientific Development Department 

Bechtel is recognized in the United States as a builder for industry 
and government, and for its many years as an international engineer, con- 
tractor and project manager of facilities for the natural resources and 
energy industries. Under single responsibility the Company provides com- 
plete service for the realization of projects, including economic and 
technical studies, engineering design, procurement, and construction or 
project management. During early years the organization participated in 
construction milestones as Hoover Dam and the San Francisco-Oakland Bay 
Bridge. More recently, activities have expanded to encompass such fields 
as nuclear engineering, refinery systems, food processing, chemical plants, 
pipelines, ore processing and other basic private industries. 

In the public sector the firm is active in water use and conservation, 
rapid transit, water desalination and the missile and space programs. 

Indicative of its concern with technological advancement, Bechtel 
has been a pioneer in the application of nuclear energy. Starting with 
construction of EBR- 1 in 1949., the first nuclear reactor in the world 
from which useful electric power was produced, the company has been a 
major participant in the conception, evaluation, design, engineering and 
construction of many nuclear projects, including fourteen generating plants. 

In recent years when the idea of combining nuclear power with a sea 
water desalination plant was proposed, Bechtel was one of the engineering 
firms thatt:ook an active interest. Bechtel had long been acquainted with 
the use of processes for desalting makeup sea water for power plants. 
The concept of large-scale nuclear power and desalting combines was a 
technically attractive prospect. Bechtel had designed and built the 
Mandalay Steam Power Station near Oxnard, California. This was the first 
power plant in the U.S. to utilize the multistage flash distillation 
process to produce the total fresh water requirement for the plant. The 
plant had 100,000 gpd capacity and was designed and built on a subcontract 
by Cleaver Brooks, Inc. 

Background on the Choice of the M.S.F. Process 

When Dr. Lionel S. Galstaun came to work for Bechtel f s Scientific 
Development Department in 1962, he further stimulated interest in the 
combined power-water concept and had a keen personal interest in the 
multistage flash process. Soon after that date Bechtel entered into a 
contract with the Office of Saline Water to make a comparative study of 
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the three most successful government pilot plants. The study had to be 
based on a normalized cost estimate so that the costs of the individual 
facilities would be comparative. Considerable criticism existed at that 
time because it was felt by many engineers that cost estimates being made 
by many different engineering concerns were not truly comparable. 

That study, published in 1963, was entitled M A Study of Large Scale 
Saline Water Conversion Plants 11 and contained the following cost of 
product water summary (in cents per 1000 gal.)# 
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This study based on the three demonstration plants then in operation fur- 
ther confirmed the belief that the MSF process held an increasing cost 
advantage as the size was increased. In fact, it was suggested that a 
350 day/yr plant availability would be justified for the multistage flash 
process instead of the normalized 330 day used in the study. 

A picture of the 36 stage San Diego Plant is shown in Exhibit 1. 
Exhibit 2 is a flow diagram for 50,000,000 gpd based on an extrapolation 
of this idea, Exhibit 3 shows a typical cutaway evaporator vessel for the 
50,000,000 gpd design concept. 



"A Study of Large Size Saline Water Conversion Plants 1 ' by the Bechtel 
Corporation made for the Office of Saline Water, U.S. Dept, of the 
Interior. R 6c D Report 72, March 1963. 

The electrodialysis plant was based on brackish water and thus eliminates 
it from our consideration when operating on sea water and above 2 million 
GPD capacity. 
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Events Leading to Crawford's Assignment 

In January 1965 a member of the Scientific Development Department 
of Bechtel attended a course in Dynamic Programming at Stanford University 
taught by Dr. Douglas Wilde. Upon reporting this experience to Dr. Galstaun, 
Dr. Galstaun suggested that it would be a good idea for Bechtel to retain 
Dr. Wilde, on a consulting basis, to develop generalized optimization 
approaches to process design problems. Dr. Galstaun suggested that process 
optimization would be of interest to the Process Engineering and Water Re- 
sources groups of the Applied Technology Section. Dr. Galstaun further 
suggested that an initial study be made to see whether new optimization 
approaches could be applied to the MSF desalination process. 

Thus, the initial design team consisted of Mr. R.L. Clark, Mr. Rufus 
Crawford, and Dr. Wilde with Dr. Galstaun supervising the study. 
Mr. Clark was a senior member of the desalination group and served as the 
direct supervisor for Rufus Crawford. Many of the decisions needed for 
Crawford to get started were made by Clark and Galstaun based on prior 
desalination feasibility studies carried out by the Bechtel Corporation. 

It was in June 1965 that Dr. Wilde was actually retained as a con- 
sultant and a start was made to optimize the MSF process utilizing a 
computer program. Soon after this date, Dr. Galstaun was assigned to the 
New York office of Bechtel and Mr. R.L. Clark was assigned to the office 
in Mexico City. The working team then became: 

Mr. James Appleton, Project Supervisor 

Mr. Rufus Crawford, Study Engineer 

Dr. Douglas Wilde, Consultant Chemical Engineer 

Mr. Moredicai Avriel, Stanford Graduate Student, Summer employee 

of Bechtel and assistant to Dr. Wilde. 
Miss Marilyn Martin, Computer Programmer 
Mr. John Chapman, Computer Programmer/Chemical Engineer 

Crawford recalled the comparative study made by Bechtel on the three 
OSW demonstration plants. Crawford said that for convenience the MSF 
process could be called the long tube horizontal (LTH) process in contrast 
to the long tube vertical (LTV) process. He said, n The processes do have 
a number of differences other than the fact that the tubes are horizontal 
or vertical." The LTV was the only other distillation plant considered 
in the comparative study. The one million gallon per day plant at 
Freeport, Texas, is of the LTV type and operating data are available. 
In this multiple effect distillation process sea water enters the first 
effect and is heated by an external source. At Freeport this is a steam 
supply from the immediately adjacent Dow Chemical plant. Steam that is 
formed in the first effect is condensed in the second effect which evap- 
orates more sea water (at a reduced pressure). This process is repeated 
in succeeding effects, the pressure changes corresponding to the temperature 
difference required to produce the flow of heat needed. See Figure 1. 
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In optimizing the MSF process, Crawford said their first inclination 
was to apply the principles of dynamic programming, a promising optimi- 
zation technique, directly to the MSF process but they found that the 
computer program required too many storage cells. They also felt a need 
to find a more efficient optimization technique capable of analyzing the 
many MSF system parameters. They knew they were going to take partial 
derivatives and this ruled out the trial and error solution analysis 
even though a cycle analysis could have been performed on that basis. 

Dr. Wilde, after studying the MSF process for about one week sug- 
gested that the Pontryagin Maximum Principle (derived in Russia) might 
be applicable to this problem and prove to be a useful tool in optimi- 
zing this system. This technique is a sophisticated application of 
variational calculus and had been investigated by Dr. Wilde at Stanford. 
However, obtaining the process variables in the proper form for the 
computer program proved to be a challenging assignment for Crawford and 
Avriel . 

"Gradually we realized, 11 Rufus said, "that our procedure could be 
divided into four stages." Rufus, himself, would develop the mathematical 
MSF model suitable for the computer. Wilde and Avriel would be responsible 
for setting up the algorithm, i.e., the mechanics of the maximum principle 
concept. The team as a whole planned to work on the adjustment procedures 
for improving the MSF model and finally, Rufus would evaluate the output 
of the computer program and the effectiveness of alternatives. The cost 
of water per 1000 gallons was to be the basic parameter to be optimized. 

The problem which was defined by Crawford in consultation with 
Galstaun and Clark was to have the following initial boundary conditions: 

1. The MSF system would have greater than 10 million gpd capacity. 

2. Sulfuric acid or polyphosphate additives were to be used for 
scale control but not to be a major optimization parameter. 

3. A single effect recycle system would be used consisting of: 

(a) heat input section 

(b) heat recovery section 

(c) heat rejection section 

These general boundary conditions were selected so as not to over- 
restrict the optimization study. Mr. Crawford had available the following 
simplified flow diagram for a flash plant using polyphosphate additive. 
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BLOW DOWN OF COOLING WATER 
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Figure 2. Simplified Diagram of MSF Process 



This cycle operated reliably at San Diego and other site locations 
throughout the world. It produced one million gallons per day when 
using the polyphosphate (HAGEVAP LP at 4 ppm) additive to prevent 
the adherence of scale (mineral deposits) to the heat transfer sur- 
faces* Then with a minimum of change they were able to produce 
1,400,000 gallons per day using a maximum brine temperature of 250°F. 
In this case, scale was controlled by the addition of (125 ppm) a 
stoichiometric amount of ^SO^ to break down the HCO^ ion present in 

sea water. The evolution of CO^ after acidifying with H^SO^ requires 

that the deareation capacity of the system be double that for a system 
which removes the dissolved gases alone. 

The range of heat input "q" for multistage flash systems varies 
from 50 to 200 Btu per pound of product water, with values near 100 
being more common in current practice. Another useful parameter, the 

gained output ratio (g.o.r.) is defined as where the 1000 is 

q 

approximately the latent heat of vaporization. For a typical one million 
gpd plant such as San Diego 



Some other conditions had to be specified in order to analyze the 
"phosphate 11 cycle specifically. The steam supplied to the heater was 
assumed to be supplied at 213°F. The brine entering the heater was 187°F 
and the other assumed conditions are shown on the simplified flow diagram 
(Figure 2). 



(g.o.r.) = 10.5 



1000 (Btu/lbm steam) 
95.2 (Btu/lbm distillate) 
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The concentration factor is defined as the ratio, 

(ppm) salts in the blowdown 
C * F * " (ppm) salts in the sea water 

for MSF cycle values are generally between 1.7 and 2 as to prevent 
scale build-up on the heat transfer surfaces , The number of stages 
may be specified. A typical value for a g.o.r. of 6 might be 25 stages 
while for a g.o.r. of 12 approximately 50 stages would be used. 

The fluid velocities found within the condenser tubes are limited 
by inlet and outlet erosion problems. Depending on the tube material 
the fluid velocity in the tubes may be taken as 4 ft/sec to 7.5 ft/sec. 
These velocities have to be set for all three sections. The velocities 
are needed to estimate the heat transfer coefficients and ultimately 
the condenser area needed for each of the plant sections. 

Crawford knew that he must derive the process analysis equations 
before Avriel and Dr. Wilde could start their optimization program. 
Some method would have to be developed to measure the effectiveness of 
this program. He was mentally planning his approach as he drove home. 
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36 -STAGE LTH DISTILLATION PLANT 



Office of Saline Water 
San Diego, California 

This 36 -stage, 1,000,000 gpd flash facility began operation as a demon- 
stration plant in late 1961. The Bechtel Corporation was selected in 
196Z by the Office of Saline Water to evaluate the plant's performance, 
and to make recommendations, and determine the cost of water. Seen 
above are the ends of the five flash vessels, the brine heater, and, in 
the background, the boiler feed deaerator and boiler stack. Since this 
picture was taken, this plant was dismantled and shipped to Guantanamo 
Bay, where it is now producing water to meet the needs of the naval base . 



EXHIBIT 1 
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As shown in the following cutaway view of a typical evaporator vessel, 
each vessel contains condensing tubes, product water troughs, and a flashing 
brine area. The condensing tubes, carrying the recycle brine, extend through 
the internal walls, thus going from stage to stage in each vessel. The flashing 
brine flows in the opposite direction from stage to stage through orifices 
in the vessel internal walls. The product water troughs for each stage are 
connected by loop seals and orifices. 




Typical Evaporator Vessel 



EXHIBIT 3 
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THE BECHTEL CORPORATION (B) 



Process Analysis and Optimization of a Desalting Plant Design 



Mr. Crawford decided it would be necessary to set up a mathematical 
model for a single stage as the first phase of an overall flow stream model. 
After the basic energy balances were determined for the single stage Rufus 
and Mordecai would extend their analysis to combined stages and finally to 
the MSF system as a whole. They expected to select variables suitable for 
the optimization scheme. 



Single Stage Model 
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Figure 3 Simplified Flow Model of the i Stage 



let X denote any property of a stream entering the stage and X that 
same property as the stream leaves. The subscripts R, D and B represent 
recycle brine, distillate and flashing brine, respectively. Crawford said 
that Dick Clark didn't like the X tilda ( X ) symbol, saying "What happens 
if you lose a "tilda"? Actually a single symbol could have been used with 
the subscript i indicating the interface between stages instead of the stage 
number itself. However, the program was already underway so the notation was 
not changed. 
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THE BECHTEL CORPORATION (B) 



Process Analysis and Optimization of a Desalting Plant Design 



Mr. Crawford decided it would be necessary to set up a mathematical 
model for a single stage as the first phase of an overall flow stream model. 
After the basic energy balances were determined for the single stage Rufus 
and Mordecai would extend their analysis to combined stages and finally to 
the MSF system as a whole. They expected to select variables suitable for 
the optimization scheme. 
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In the Simplified Flow Model of the i Stage, Avriel and Crawford 
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Figure 3 Simplified Flow Model of the i Stage 



let X denote any property of a stream entering the stage and X that 
same property as the stream leaves. The subscripts R, D and B represent 
recycle brine, distillate and flashing brine, respectively. Crawford saxd 
that Dick Clark didn't like the X tilda ( X ) symbol, saying "What happens 
if you lose a "tilda"? Actually a single symbol could have been used with 
the subscript i indicating the interface between stages instead of the stage 
number itself. However, the program was already underway so the notatxon was 
not changed. 
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Since each stage is assumed to be at a constant pressure, (flashing 
takes place at the stage interface as the brine is throttled through an 
orifice) the relationships between stream temperatures may be deduced from 
an enthalpy balance. Thus the enthalpy lost over the stage by the brine 
stream must equal the enthalpy gained by the distillate and recycle brine 
streams. Crawford made the decision as to which brine stream heat capacity, 
density, thermal conductivity, viscosity and boiling point elevation data 
would be used in the calculations. He mentioned that it was accepted practice 
by now to use the data from the "Saline Water Conversion Engineering Data Book 11 
compiled by the ML W. Kellogg Company under contract with the Office of Saline 
Water of the U. S. Department of the Interior. Crawford did not check these 
data directly but the known thermodynamic properties of fresh water and 
extensive data on NaCl solutions serve as a rough check. These data may be 
found in Exhibit 4 originating from Oak Ridge National Laboratory., Crawford 
pointed out that, "The method of solution is independent of the data fed 
and a higher degree of accuracy can be achieved as more reliable data become 
available." For purposes of checking the single stage heat balance equations 
the variations in brine properties with concentrations was unnecessary. 

Combined Stage Model 

After completing and checking the single stage analysis, Crawford and 
Avriel expanded their concept to include the entire heat recovery section 
(see Exhibit 5) and with Miss Martin's assistance programmed the equations. 
They analyzed a 33 stage operation so that a check might be made with an 
earlier report, Bechtel No. 4427 which was issued by the Office of Saline 
Water as R, 6c D, Report No. 116 and may be purchased from the Superintendent 
of Documents, U, S, Department of Commerce. This report used a different 
computational sequence and was somewhat lengthier. This comparison saved 
them the job of making a hand calculation to check the computer solution. 
The results of the analysis showed that the new relations used by the computer 
appeared to be consistent with those used in 4427 but a final judgment was 
not made at this time since the heat transfer coefficients were slightly 
different. 

The next step was to write a heat balance on a single stage of the heat 
rejection section. The rejection stages were then combined. The brine 
heater was then included assuming inputs from a steam source of + 8° F 
above the maximum allowable brine temperature. Crawford said that the program 
will accept any A T but that Mordecai did a detailed analysis of this using 
geometric programming and concluded that 8° F was very near the optimum for 
this specific application. The recovery, rejection and brine heater sections 
were then combined to form an integrated system to which boundary conditions 
could be applied, 

The Avriel letter (Exhibit 5) records the heat and material balance 
relations of the MSF study. Crawford said, "I would have liked this material 
to have been better organized for the record, but I agree completely with the 
conclusions." In fact, he felt he was a part of those conclusions. 
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Part of the study had been to derive equations for the heat transfer 
coefficients. Avriel used the standard textbook equations such as the 
Dittus^-Boelter equation. These had been used by Bechtel prior to the study 
and had been shown to reproduce with reasonable accuracy heat transfer data 
taken from the San Diego and other facilities. They were eager to get the 
results of the first computer run since they were not 1007 o sure that they 
hadn't made some mistake in their derivations and wanted to check out the 
program. They were using report No. 4427 but they could have used actual 
operating data such as that of the San Diego Plant (Exhibit 6) . 

The cycle analysis was recognized as an application of the first law 
of thermodynamics and since this is true by definition, what goes in should 
balance what comes out. However, there was an error somewhere because the 
computer program repeatedly gave an imbalance. This showed up as a predicted 
sea temperature of 61° F instead of the 60° F temperature originally postulated. 
The cycle balance was accomplished by writing equations for the first stage 
(stage nearest the heater); then the outputs from this stage were fed to the 
second stage until a balance was achieved then fed to the third, etc. The 
final values obtained in this manner should have checked with an overall 
balance for the system. However, they did not and the error indicated that 
more heat was going in than was coming out. In order to place a high degree 
of confidence in the program the inconsistency had to be resolved. Mr. Crawford 
saw three possible sources of error. First, his hand calculation might have 
been in error. Or, the basic equations may not have been set up correctly. 
Finally, some error might have occurred in the computer program. 

Mr. Crawford started the tedious task of running down this inconsistency. 
He rechecked his calculations using a slide rule but soon realized that the 
significant figures given by the slide rule would yield a greater variation 
than the discrepancy for which he was looking. 

Next he checked his equations and could find no error. He questioned 
the accuracy inherent in the method of calculation and concluded that this 
could account for only 1/4° F in the sea temperature which was the final 
output number. He then started to check the computations stage by stage on 
an electric calculator, using eleven numbers on the readout register so that 
the question of significant figures would not arise again. 

In order to make the length of this task reasonable and because there 
was a question in Crawford 1 s mind as to the influence of the empirical 
equation for the specific heats, Crawford assumed all specific heats to be 

B tu 

1.0 — — and reran the computer program with this constraint. By the time 

Crawford got through the heat recovery section the numbers were still checking 
the computer program. However, in going through the heat rejection stages 
the error reappeared so he concluded the error must be in that part of the 
program. 



B-3 



ECL 70 



Avriel and Wilde were also searching for the error and they had about 
concluded that the problem was overspecif ied and had gone to Dr. Galstaun to 
discuss this possibility with him. Almost at this exact time John Chapman 
called to say the program had just run. 

It appeared that the error had been in the equation for computing the 
constant A^. It seems the equation had been written in the program as 

A RN = "AVi instead of Arn " «l W w?pPi 



where W RN is the circulating cooling sea water in the heat rejection section 
and the computer program was calling for W the circulating water in the 
heat recovery section. 



Student Assignment 



(1) Set up the equations for heat and mass balance for the "i 11 stage 
in the heat recovery section. (as Crawford did) 

(2) Set up the equations for heat and mass balance for the M i th " stage 
in the heat rejection section. (as Crawford did) 

(3) How should Crawford now proceed? 
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Exhibit 4 



Physical Properties 

Physical properties will be calculated from the following relations (Oak 
Ridge) 

1. Boiling Pt. Elevation (Low Concentration) 

Concentration range: .02 - .12 wt. fraction, solids 
Temperature range: 77 - 392°F 

A(°F) = A1+A2(C)+A3(T)+A4(CT>A5(C) 2 +A6(T) 2 +A7(C)(T) 2 +A8(C) 2 (T) 

+A9(C) 3 +A10(T) 3 +A11(C)(T) 3 +A12(C) 3 (T)+A13(C) 2 (T) 2 +A14(C) 2 (T) 3 



where 



C = wt. fraction, solids 

T = temp. , °F 

Al = 0.1014380 E-l 
A2 = 0.1021815 E+2 
A3 = -.6386588 E-4 
A4 = 0.6118005 E-l 
A5 = 0.8809554 E+l 
A6 = 0.2386759 E-6 
A7 = 0.2214495 E-4 



A8 = 0.1714722 
A9 - 0.7173776 E+2 
A10 = 0.1170617 E-9 
All = 0.1042795 E-6 
A12 = 0.3620461 
A13 = 0.7823281 E-4 
A14 = -.5218751 E-6 



2 . Boiling Pt. Elevation (High Concentration) 

Concentration range: .12 - .28 wt. fraction, solids 

A(° F) = A1+A2 (C)+A3 (T)+A4 (C) (T)+A5 (C ) 2 +A6 (T) 2 +A7 (C) (T) 2 +A8 (C ) 2 (T) 

+A9(C) 3 +A10(T) 3 +A11(C)(T) 3 +A12(C) 3 (T)+A13(C) 2 (T) 2 +A14(C) 2 (T) 3 



where 



+A15(C) 3 (T) 2 



C = wt. fraction, solids 
T = temp., °F 



Al 




-.1441403 


E+l 
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= 0.4579313 E+3 
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E+2 
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= 0.4010406 E-7 
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0.1081565 


E-l 


All 


= -.4393590 E-6 


A4 




-.1253869 




A12 


= -.1505255 E+l 


A5 




-.1679212 


E+3 


A13 


= -.3665256 E-2 


A6 




-.4486056 


E-4 


A 14 


= 0.1280709 E-5 


A7 




0.7590168 


E-3 


A 15 


= 0.5852867 E-2 


A8 




0.1200749 


E+l 
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3. Specific Heat : Btu/(lb )(°F) 

Data available for 40 to 250°F and 0 to 26% salt. The data 
was extrapolated to 400°F. 

C p - 1.0057733-.13280442X10~ 3 T-1.5344907C+. 39090715 x 10 _2 TC 

+,65092605X10' 6 T 2 +2.4955446C 2 - . 11156771xl0" 1 TC 2 

-.88941836xl0" 5 T 2 C+.3598702xl0" 4 T 2 C 2 

where C - wt. fraction of salt 
T - °F 

3 

4. Density : lb /ft 
m 

Data available from 40 to 212°F and 0 to 26% salt. The data 

has been extrapolated up to 400°F using data of pure water as a guide, 

p - 62.717753-.32152986xl0" 2 T +.44315006X10 2 C- . 11647394X10" 1 CT 

-.48932777X10" 4 T 2 +.16449945X10 2 C 2 +.13761984X10" 1 C 2 T 

+.34231326X10" 5 CT 2 -.92275301xl0" 4 C 2 T 2 

C - wt. fraction of salt 
T - °F 

5. Thermal Conductivity: Btu/(hr)-(Ft)(°F) 

Data available from 40 to 250°F and 0 to 26% salt. 

Data has been extrapolated to 400°F using purewater data as a guide, 
k - -.0655C + .29556783 + .8794748X10" 3 T 

-.23026302X10" 5 T 2 + . 1609445X10" 8 T 3 
c - salt cone. (wt. fraction), T-°F 

6. Viscosity : LB /(hr)-(Ft) 

Data available from 40 to 250°F and 0 to 26% salt. The data was 
extrapolated to 400°F 

\i - 22. 216306-14. 161033C-16. 460588 (log 1Q T ) +16.430566 (log 1() T ) C 
+284. 84684C 2 +3. 0949585 (log 1() T ) 2 -4.205263C (log 10 T ) 2 
-219.40291 (log 1() T )C 2 +42.287729C 2 (log 1() T ) 2 

C - salt cone. (wt. fraction), T-°F 

[i - (lb /hr-F T ) 

m l 

M. Avriel 
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BECHTEL CORPORATION 

INTER-OFFICE MEMORANDUM 

Date July 19, 1965 

To M. G. Martin From M. Avriel 

Dept. Scientific Development 

Copies to L. S. Galstaun At 301 Mission 

D. J. Wilde 

R. L. Clark Subject: 

R. W. Crawford 




Programming Heat & Mass Balance Relations of 
MSF Optimization Study 

We are now ready to program and test the heat and mass balance 
relations for the optimization study. These relations determine 
the physical conditions at the outlet of each stage as a function 
of the physical conditions at the inlet to the same stage. Knowing 
the terminal conditions at one end of the flashing stages enables 
one to calculate recursively the conditions at each stage. 
Let 

Hr^ = enthalpy of recycle brine, entering stage i, (Btu) 
Hb^ = enthalpy of flashing brine, entering stage i, (Btu) 
H n . = enthalpy of distillate, entering stage i, (Btu) 
H Pi = enthalpy of distillate produced in stage i, (Btu) 
Tr^ = temperature of recycle brine entering stage i, (°F) 
Tb-£ = temperature of flashing brine entering stage i, (°F) 
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T 

Di = temperature of distillate entering stage i, (°F) 

= flow rate of flashing brine entering stage i, (lb/hr) 
= flow rate of distillate entering stage i, (lb/hr) 
= flow rate of recycle brine in heat recovery section 

C p = specific heat (Btu/lb°F) 

X = latent heat of vaporization (Btu/lb) 

A = boiling point elevation ( °F) 

2 

U = overall heat transfer coefficient (Btu/hr ft °F) 
n^ = number of tubes in stage i 

= 0.D, of tubes in stage i (ft) 

= tube length in stage i (ft) 

Also let 

denote a variable leaving stage i, or a physical property corresponding 
to the conditions of a stream leaving stage i For example Cpg^ ^ s tne 
specific heat of the flashing brine at temp and concentration as it leaves 
stage i. 

X^ denote an average variable throughout a stage unless specified otherwise. 

The following basic relations exist at each stage of the recovery section 
Overall heat balance on stage i 

~m 

( 1 ) . + H . + H . = H + H + H 

Di Bi Ri Di Bi Ri 

Heat transfer at stage i. 



T Ri L 



(2) 

T 



/ \ Sr dT R " \\ n i « °i ( " D i - T R } dX 



Ri 
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where it is assumed that the temperature of the distillate is constant 
across the stage and it is equal to the equilibrium temp. 
Heat balance on flashing brine: 

(3 > =81 - «Bi + «Pi + V W Bi ""El > 

Define the following average quantities : 
Average specific heat of distillate at stage i: 

(4 > - : c pdi + Sdi 

PDi 2 
Average specific heat of recycle brine at stage i: 

(5 > jf :.. J Sri + Sri 

PRi 2 
Average specific heat of flashing brine at stage i: 

< 6 > - ; ' c PHi + ~pbi 

PBi 2 

Note that (4) and (5) are functions of inlet and outlet temperatures, while 

(6) is a function of inlet and outlet temperatures and concentrations. 

Weighted average overall heat transfer coefficient at stage i: 

(7) _ J U. ( ? M -T R . ) f \ ( ? D . -~ R , ) 

-i mm mm mm 

(T-T) + (T-T) 
V ^Di R ' k X Di Ri ; 

Average boiling point elevation at stage i: 

mm 

(8) _. A. + A. 
A i = 1 2 1 

Eq. (8) is a function of both temperature and salt concentration 
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Let 
(9) 



a = 

i 



P. = 

1 



/ C: . A. 
i/ PBl i 



exp 



- 1 



w r Sri 



Then, combining eqs . (1) through (9) one obtains 

f 



a.W„C m .P.l . + J a.W„. (C^ . - C, 



i R PRi i Ri 



+ a.W_. 
| l Di 



PBi PDi 



) + 



PBi 



Di 



a.W C _ .p. + J a.W_. (C_ B 
L R PRi L l l Di PBl 



C PDi> + [ W 



-(a. + l)W R ] 



J PBi 



The above equation does not take into account conversion of pumping work into 
sensible heat. 



P. S a .^. (C^ . - C .) + [w^. - (a. + 1) wl 
i J i Di PBi PDi' [ Di l Rj 



C PBi (T Di " ? Ri> 



= T Ri + 



a.WC^.P. + a W . (C . - C TO .) + IV. - (a. + 1) W 1 C: . 
i R PRi i i Di PBl PDi [ Di 1 R J PBl 

(W„ - W^.) axL.P. (T^. - tV) 
R Dr i R PRi i Di Ri 



= w^. + — r 



a. A. i a.W^C^.p. + aW (C TC . 

i Di PBi 



v i~i | i"R~PRi 1 
Note that by defining 
(13) 



. - C^.) + fw^. - (a. + 1) W_l 
l PDi' [ Dl l Rj PBl 



(14) 



A„ = a.W C B. 

Ri l R PRl l 

B„ = a.W„. (C TO . - £.) + - (a. + 1) W 1 C__ . 

Ri l Dl PBl PDl I Di l Rj PBl 



eqs. (10) to (12) become respectively: 



(15) 



T . + B„ T . 
Ri Ri Ri Di 



Di 



. Ri Ri 



(16) 



P.B„ . (T^. - T .) 
i Ri Di Ri 



T . = T . + 
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(17) w - ^. 4 

Di Di 



(W R * W Di ) 



a. A. 



A„. (T - T.) 
Ri v Di Ri' 



A Ri + B Ri 



Eqs. (15), (16) and (17) are the basic equations for stagewise 

calculations. However, they cannot be evaluated in a straightforward 

manner since A and B both contain average values of physical properties. 
R R 

Therefore, a converging iterative procedure is required to calculate 
eqs . (15) to (17) simultaneously . 

Let the stages of the heat recovery section be numbered 1, 2, N, 
starting at the brine heater. The heat rejection stages are numbered N+l, 
N+2, ... For the heat rejection section the following assumptions are made: 

1. Cooling water and make-up water may flow through separate 
tube bundles with different temperatures. 

2. The distillate leaving stage N may be withdrawn from the 
system. 

Let 

T t/ xt i • = temp, of make-up water entering stage N+i ](°F) 
MjN-ri 

T ■ - temp, of cooling water entering stage N+i (°F). 

Cj , JN"rl 

W = flow rate of make-up water (lb/hr) 
M 

W r = flow rate of cooling water (lb/hr) 

The meaning of all other variables appearing in the following relations is 
consistent with previous discussion or should be clear from the subscripts 
used . 
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In addition we define: 
(18) 

P M,N+i 



(19) 



where 
(20) 



1 - exp 



C,N+i 



1 - exp 



11 irD L n 

M,N+i * M,N+i N+i M,N+i 



W M C PM,N+i 



U c,N+i rtD c,N+i L N+i. n c,N+i 



W C KU-i 
c pc , N+L 



U, 



M,N+i 



U M,N+i (T D,N+i " T M,N+i ) + U M,N+i (T D,N+i " T M,N+i } 
(T D,N+i " T M,N+i ) + tT D,N+i " T M,N+i ) 



and U has the same form as (20), except for replacing M by C everywhere 

c,N+i 

at the right hand side- 

If the make-up and cooling streams have identical temperatures and flow 
conditions then eqs. (18)-(20) are not necessary and eqs. (15), (16) and 
(17) are used with the following changes: 



(21) 
(22) 



W R — > W RN - W M + W c 



> T R,N+i T M,N+i T c,N+i 



For the general case and for no withdrawal of distillate at stage N: 



(23) 



D,N+i 



AT +AT +BT 

A M Si, N+i c c,N+i R D,N+i 



A M + A c + B R 
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where B Rn - a^. (C pB . - C pD .) + 



[ W Di " ( «i + X) \] 



PBi 



(24) 
(25) 

(26) 
(27) 



A ■ a W C 3 
M i M PM,N+i H M,N+i 



a.W C_ , . 3 , . 
i c Pc,N+i c,N+i 



If make-up and cooling water streams are identical: 
A,. + A = A, 



M c 



RN 



and 



T = T = T 

M,N+i c,N+i R,N+i 



Similarly, 



(28) 



(29) 



M,N+i 



c,N+i 



(30) 



W. 



D,N+i 



T . + 
M,N+i 



T . + 
c,N+i 



W D,Iffi + 



P M,N+i A c (T c,N+i " T M > W-i ) + P M,N+i B RN ^ T D,N+i " T M,N4-i ) 

A M + A + B_ M 
M c RN 

*+> *» 

P c,N+i A M (T M,N+i " T c,N+i ) + ^c.wAn (T D,N+i " T c,N+i ) 



(W R " ^D,N+i^ 



A H + A c + B RN 



(A M + V T Di " VM.N+i " A c T c,N+i 



A + A + B 
M c RN 



Here eqs. (23), (28), (29) and (30) are the basic equations, corresponding 
to (15), (16), and (17) above. 

If the distillate is withdrawn after leaving stage N, the following 
relations are valid in the heat rejection section: 



(31) 
(32) 



W D,Nfi * ° 



B,N+i W B,N+i " W D,N+i 
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Let 



(33) B BN + a N + i W D,N + i ( V N+ i - C PD,N +i ) + 



[ W D,N + i " ^ 



N + i + X) W B,N + i 



PB,N+i 



Eqs. (23), (28), (29) and (30) are used then with B replaced by B 

RN BN e 
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Heat Transfer Calculations 

Heat transfer coefficients will be calculated at conditions corresponding 
to both ends of a stage and an average will be taken. 

Assumptions : 

1, Constant vapor temperature, T^^ across stage i. 

2, The condensing film temperature has a greater effect upon 
the tube wall temperature than the bulk recycle temperature. 

3, Temp, drop across tube wall negligible for purposes of wall 
temp . calculations . 

The following wall and film temperatures will be used: 
1, Wall temperature 



T . . - 0.7 T n . + 0.3 T_ . 
Wi Di Ri 



and 



T. = 0.7T+O.3T. 
wi Di Ri 



Condensing film temperature: 



3 1 ^ 

T — T + — T 

Dfi 4 wi 4 Di 



3 ~" 5 ~« 

t = — T + — T 

x Df i 8 Ri 8 Di 
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Df i 



3 5 

— T + — T 

8 A Ri ^ 8 Di 
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The following film coefficient correlations will be used: ( "Transport,. 
Phenomena" Bird, Stewart, Lightfoot, Wiley p. 399) 

1. Inside the Tube 



h.,. = 0„026 
Ri 



k . 
Ri 

(ID)j 



Ri 



l"l 




1/3 


1^1 


0.14 













where k Ri = thermal conductivity of recycle brine, Btu/(hr)(sq, ft)(°F)/ft, 

viscosity of recycle brine, lb/(hr)(ft) 



Ri 



Wi 



Ri 



U Ri 
<I.D.) 4 



viscosity of recycle brine evaluated at the wall temp . , lb/ (hr) (f t) 

3 

density of recycle brine, lb/ft 
velocity of recycle brine, ft/hr 
inside diameter at stage i, ft 



All physical properties are evaluated at T or T 

Ri Wi 



S imilar ly, 



h B . = 0.026 
Ri 



Ri 
(I.D.). 



Ri 



C U 
PRi Ri 



1/3 



"ill 



^Ri 



0.14 



where the variables have the same meanings as before, but calculated at temps. 



corresponding to TL . or T„.. 

Ri Wi 
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2 . Condensing film coefficient 

The following correlation will be used: ( "Unit Operation 11 , Brown, p. 449) 
( k DfJ IpDfJ g V 



h Di 



0.725 



N. D. u_. (T_. - T ) 



where 

X i = latent heat of vaporization at T^ 

All other physical properties are calculated at temp. 

N. = no, of tubes in a vertical column of tubes. 

L 



Similarly, 



h D . = 0.725 



r .3 , ,2 ~ ,1/4 

( k Dfi) Upfil g Y 



N i D i %fi (T Di - V 



Overall heat transfer coefficient, based on outside diameter: 

1 



U. = 

1 



+ [ D i - (TD) i ] "i 

h Di h Ri (ID) i ( D i +(ID) i ] ( k m) 



+ FF 



where k « thermal conductivity of tube metal, Btu/(hr)(sq. ft)(°F)/ft 
m 

FF - fouling factor, (hr) (sq . f t) (°F)/BTU 
Similarly, 



U i * 



_ + K - (iD >iK D i> 

VV ID) i + [ D i +(ID) i] ( k m) 



+ FF 



FF was linear with temperature, °F, varying from .0005 at 85°F to 



0.001 



ft ,hr,°F 
Btu 



at 250°F. 
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BECHTEL CORPORATION (C) 
Process Analysis and Optimization of a Desalting Plant Design 



After studying the multistage flash process Dr. Wilde suggested 
that the "Pontryagin Maximum Principle 1 ' might be a useful optimization 
tool. This method is a sophisticated application of variational calculus 
and had been examined by Dr 0 Wilde at Stanford . Crawford and Avriel 
found it necessary, as in most optimization processes, to develop explicit 
relationships for the process variables from the cycle analysis so that 
partial derivatives could be conveniently determined, Crawford summarized 
(Exhibit 7) his energy and material balance equations for the computer 
program. 

Initially Rufus decided to base his computations on cylindrical outer 
shell evaporators and worked with this geometry for about one month . 
Since previous experience at San Diego indicated that a rectangular shape 
allowed a larger volume of brine flow, Clark and Galstaun suggested that 
Crawford switch geometries. The original rationale for choosing a 
cylindrical cross-section was based on the loading of the shell when 
under vacuum. However, the magnitude of the loads actually encountered 
were such that no serious stress problems developed in the San Deigo 
rectangular shells. Some question existed as to whether the cylindrical 
or rectangular shape was the more economical evaporator configuration . 

The cost estimates for this program were based on steel pressure 
vessels . Earlier, Bechtel had made a study comparing coated concrete 
and steel shells. Although the cost estimates indicated a competitive- 
ness for concrete, limited experiments indicated that the surface coatings 
would not stand up under exposure to hot brine water . It was assumed in 
the study that if concrete were used for the pressure vessels, the slabs 
would be cast at the construction site. Crawford said, "At the present 
time the absence of data on the concrete application makes it difficult 
to predict the long-term reliability of these evaporator vessels," 

The tube material and its thermodynamic properties were inputs to 
the computer program. Aluminum-brass was used since this provided a 
check on the 4427 study. However, based on power plant experience the 
general feeling is that if they are going to design for a 30 year life 
that 70-30 Cu-Ni will have to be used „ Crawford said, "Cu-Ni is probably 
the best material you can afford to use." The tube diameter is an input 
to the computer program and 5/8, 3/4, 7/8 and 1 inch are standard, and 
currently the maximum mill tube length is 60 feet. The 3/4 inch diameter 
(O.D.) is most common and is usually assumed for the first try. 
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The question of deareation arose repeatedly but the decision was 
made to include it only as a fixed cost item„ Crawford though to himself, 
"This would be an optimization study in itself . 11 In estimating the cost 
of deareation it was included as an integral portion of the heat rejection 
section , 

For scale control Crawford was told by Galstaun to design so that 
either polyphosphate or sulfuric acid could be used and these would set 
the temperature limitations of 195° and ?50°F, respectively . Crawford 
says that if anyone comes up with a new method of scale control that 
would remove these temperature limitations, new cost data could be fed 
into the computer program. As for the inherent drop in performance with 
time, this was taken care of as a fouling factor. The fouling factor is 
assumed to be a function of the temperature and the condenser tube ma- 
terial as described in part (B)„ 

On the question of drives, Crawford said he made the decision to 
limit the study to electric drives since cost data were more readily 
available than for other systems and this seemed to have the approval 
of Galstaun. Steam drives were not considered even though they were 
known to have been used on the main pumps of the San Diego plant, 
Crawford decided to wait until data on the steam drives became available 
before incorporating it into their optimization model* 

Crawford discussed the philosophy of cost estimation. He noted 
that the n broad brush" treatment often gave answers within 5 to 10% of 
detailed cost analysis and in much less time. Therefore since the 
refined equations are proprietary ,Ruf us offered some of the broad brush 
guide lines for establishing costs, 

Crawford noted that steam costs average from 20 to 30 cents per 
million Btu at 250°F and would be more expensive at higher temperatures . 
However, to lower the steam temperature at the brine heater more surface 
area is required , 

In sizing heaters and other extensive equipment he said it's pos- 
sible to roughly estimate material costs on a lineal or square foot basis. 
"For example, you can figure something under $1.00/ft^ Q f heat transfer 
surface to which must be added the cost of the tubes $2.80/ft^ Cu-Ni or 
$2,30 Al-Brass," said Crawford „ A large brine heater may then run 

$4;50/ft^ if the foundations, piping and removal of non-condensibles 
are included. The vessels which contain the tubes might run from $20 
to $30 per ft^ of evaporator bottom surface area including auxiliaries. 
A typical height of these vessels is 8 feet, 40 to 60 feet long and 
50 feet wide „ 
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Crawford, off the top of his head, set down in this fashion many 
direct and indirect operating costs. A typical breakdown (Exhibit 8) 
is reproduced from Bechtel's report to the Office of Saline Water 
(OSW No. 116). In a similar manner for any given set of boundary con- 
ditions (and with greater accuracy by computer input), Rufus was able 
to obtain a final cost for producing water given as cents per thousand 
gallons . 
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Assignment 

1„ Outline how you would attack the problem of process optimization, 

2o Crawford stated that there are a number of approaches to cycle analysis 3 

(a) Assume equal temperature drop per stage. 

(b) Assume equal condenser area per stage . 

(c) Assume equal distilled water production per stage „ 

(But that it doesn't make much difference which one you use) 
How do you think Crawford came to this conclusion? 
How else might you get started? 

3 0 From the data in Exhibit 8 derive cost relations for use in the 

computer program. If these data are not sufficient, what additional 
data did Crawford search for? 
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BECHTEL CORPORATION 

INTER-OFFICE MEMORANDUM 



To M. G. Martin 



Copies to L. S. Galstaun 
R. L. Clark 
M. Avriel 



Date November 3, 1965 

From R. W. Crawford 

Dept. Scientific Development 

At 301 Mission 

Subject: Final Relations for Heat 
and Material Balances for 
Optimization Project 
Program. 



I have found it necessary to include one additional input for 
solution of the MSF overall heat and material balances. The memo- 
randum that was sent to you dated November 2, 1965 should be 
destroyed. The following procedure should be programmed. 

1. Inputs to program: 

a. Plant Performance Ratio. 

b. Brine Temperature at Heat Outlet (Tgi) 

c. Product Flow Rate and Temperature (Wp, Tp) . 

d. Concentration of Recycle Brine (Cr) 

e. Temperature Difference Between Exiting Makeup and Exiting 
Flashing Brine (ZITm) . 

f . Temperature Difference between Exiting Coolant and Exiting 
Flashing Brine (^T c ) . 

g. Concentration and Temperature of Sea (C g , T g ) . 

2. Using values for T B? i, Tp, Wp and C R , calculate C^k, T b ^ k and 
W R as described in an earlier memorandum on recycle rate calcula- 
tions. 

3. With the plant production (Wp) , calculate the total heat duty at 
the brine heater vs ing the performance ratio given for the plant. 

q = Wp/PR 

where: PR = Performance Ratio (lb. Product/1000 Btu) 

4. Calculate the brine temperature entering the brine heater (T R ^) 
as follows: 

a. With T B5 i and C R , calculate the heat capacity of the brine, Cp. 

b. Solve for T R i (exiting recycle temperature from stage 1) as 
below. 
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t:. = t 



R ;i b,i w R c F 

c„ With C and (T +T )/ 2 (the average brine temperature 
R d 3 1 R, 1/ 

in the brine heater) calculate a new heat capacity of the 

brine (Cp) « 

d, Repeat Steps b and c until the difference between values 

of T - on two successive iterations ^p„005 o F < , 
R 3 1 

Calculate the blowdown flow c 

Wp 



6. Calculate the makeup flow. 



w M = wp + w BL 

7. Calculate the enthalpy of the blowdown initial recycle as: 

W„ Tj ,h^ + W„h w -r 

h = h = MM ' J 
BL R,N W BL + W R 

8. Substitute h . into the heat recovery section heat balance to 

R 3 N 

solve the total energy leaving the heat recovery section with the 
distillate and brine streams . 

(VjVj + W B,J h B,j) = W R ( h R,N + h B,l " Yl) 

9. With AT and T __, calculate T_ and h^ T . 

L DjK C , J C , J 

10. Substitute values obtained into heat rejection section heat balance 
obtain coolant flow rate, 

w : - h R.N (W BL + V - (W D,J h D, J ^B,J h B,J> ~_ Vgg 

c h s " h c,j 

11. Start heat recovery section stagewise calculations. Terminate th« 
section when the recycle temperature for a stage i equals T R ^. 

Test necessary for each stage calculated . 
12 - Calculate heat rejection section. 
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It has been suggested that a possible alternative to specifying an 
outlet coolant temperature is determination of a coolant flow based upon 
the size of the recycle flow. That is: 



This alternative should be programmed . The input to the program would 
not require AT , but would instead require the percentage factor to be 
used for calculation of the coolant flow. The following changes will 
be necessary in the basic calculation procedure: 

9. (Alternate) With the percentage factor, calculate the coolant flow 
rate knowing the recycle flow. 

10. (Alternate) Substitute values obtained into heat rejection section 
heat balance to obtain outlet coolant enthalpy. 



h 



w M h s - wphpf (w D)J h D>J 




L R,N (W BL + V + hS 



C, J 



w 



c 




R.W. CRAWFORD 

Scientific Development Department 



RWC : pg 



c-1 
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OPERATING COST 



Capacity ; 

Investment ; 
Load Factor; 



Power 134 mw, net 



water 50xlO c 
$60,910,000 
90 . 4% 



GPD 



A. Direct Operating Cost 

1. Fuel 

2. Chemicals 

Sulfuric Acid 

Chlorine 

Others 



B. 



Units 



10 6 Btu 



Tons 
lbs 



3. 



4. 



Supplies and Maint. Material 
(0.0015% of Invest. /Day) 

Subtotal 



$/Unit 

0.20 

20.00 
0.0635 



Labor 

Operating, Power Men 27.96 
Water Men 27.96 
Maintenance (0.0015% of Invest/Day) 
Total Direct Labor 
Payroll Extras, 15% of Direct Labor 
Total Labor Payroll 
TOTAL DIRECT OPERATING COST 

Indirect Operating Cost 

5. General and Admin. Overhead (30% of Payroll) 

6. Amortization (0.0224% of Capital) 

7. Taxes and Insurance (0.006% of Capital) 
TOTAL INDIRECT OPERATING COST 



Units/D 

51,240. 

65.7 
250. 



15. 
12. 



$/0p. Day 

10,248 

1,314 
16 
1 

914 



12,493 



419 
336 
914 



1,669 
250 
1,919 
14,412 



576 
13,600 
3.650 
17.826 



C. Total Direct and Indirect Operatint Cost 

8. Interest on Working Capital (0.725% of "C") 
TOTAL OPERATING COST 

D. Electric Power Credit, 134 mw @ $3.1l/mwh 



32,238 
234 
32,472 

10,000 



E. Net Cost of Water 

Net Cost of Water $ per Thousand Gallons 



22,472 

0.449 



*Reproduced from Office of Saline Water Report No. 116 (Oct., 1964) 
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ESTIMATE SUMMARY 



ITEM 



I Materials 
Subcontracts 
Labor 

Total Direct Cost 

II Contractor's Field Office Services 

III Startup Cost 

IV Total Field Cost 

V Contractor's Home Office 
Services, Incl. Engineering, 
Procurement and Management 

VI Land for Plant and Canal 

VII Contingency 

VIII Total Construction Cost 

IX Interest During Construction 

X Total Plant Investment 

XI Total Plant Investment (Both Plants) 
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Thousands of Dollars 



PLANT NO. 1 



Water 


Pover 


Notes 


20,838 


10,495 


CD 


3,620 


6,190 


(2) 


1,365 


1,795 


(3) 


25,823 


18,480 


(4) 


765 


920 


(5) 


150 


150 


(6) 


26,738 


19,550 




3,600 


3,000 


(7) 


60 


20 


(8) 


3,000 


1,500 


(9) 


33,398 


24,070 




2,002 


1,440 


(10) 


35,400 


25,510 


(ID 



60,910 
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WATER PLANT 
ESTIMATE IN $M 





Quantity 


Materials 


Subc'ts. 


Labor 


Total 


Evaporators 












Concrete Structure 


36 Stages 




$ 2,520 






Condenser Bundles 


7 ,050,000 SF 


$ 14 , 800 




$ 286 




Internal Piping & Baffles 


432 Sets 


1,540 




14 1 




Insulation, Lining & Paint 






/IP 

415 






Additions for Heat Rejection 


2 Stages 


100 


50 


13 




Subtotals 




$16,440 


$ 2,985 


$ 440 


$19,865 


Sulfuric Acid Tank 


10,000 bbl 




20 






CCL Release Tank 


50' x 200' 




75 






Chlorine Tank and 












Accessor ies 




/ u 




0 


117 

XX/ 


Brine Heaters (3) 


75,600 SF 


302 




2 




Product Coolers (3) 


34,800 SF 


139 




1 




Condensers (3) (Ejector ) 




40 




1 


. 485 


Recycle Brine Pump 


15,000 HP 


1,050 




10 




Other Pumps (14) 


7,186 HP 


620 




15 


1,695 


Ejectors 


3 Stages 


15 




1 


16 


Air Compressors and 












Accessories (2) 


1,200 SCFM 


50 




4 


54 


TOTAL MAJOR EQUIPMENT 




$18,656 


$ 3,080 


$ 476 


$22,412 


Instruments 




$ 165 




$ 3? 


$ iy / 


Piping 6c Valves (External 




1,000 


O A 


ion 
loU 


1 O 1 A 


Intake Structure 




1 f A 

150 




125 




intaKe otrucLUie iLCjuipment: 












and Controls 




200 




40 




Outlet Canal 






60 






E jector Saructures 




7 




2 




Pipe Support Structures 




25 




25 


634 


Insulation 






135 




135 


Electricals 




400 




110 


510 


Founda t ions ( exc 1 „Evap .Cells) 




75 




100 




Piling (Under Pump Bays) 






10 




185 


Sewers and Paving 




10 


15 


15 


40 


Painting (excL Evap > ) 






20 




20 


TOTAL - OTHER MATERIALS 




$ 2,032 


$ 270 


$ 629 


$ 2,931 


PLANT TOTALS 




$20,688 


$3,350 


$1,105 


$25,143 


Add Prorated Portion of Common 










Facilities 




150 


270 


260 


680 


TOTAL DIRECT COST 




$20,838 


$3,620 


$1,365 


$25,823 



